Possible hadronic molecule structure of 
the Y(3940) and Y(4140) 



Tanja Branz' Thomas Gutsche and Valery E. Lyubovitskij^ 

Institutfiir Theoretische Physik, Universitdt Tubingen, 
Auf der Morgenstelle 14, D— 72076 Tubingen, Germany 

Abstract. In the present article we report on evidence for hadronic molecule structures in the charmonium sector. In particular 
we discuss the K (3940) and the recently observed 7(4140) as heavy hadron molecule states with quantum numbers J — 0+^ 

or 2++. The 7(3940) state is considered to be a superposition of D*+D*^ and D*"/)*", while the 7(4140) is a bound 
state of £)*+ and D*^ mesons. We give predictions for both the strong 7(3940) — > J/xj/a, 7(4140) J/wi' and radiative 
7(3940)/7(4140) yy decay widths in a phenomenological Lagrangian approach. The results for the strong hidden charm 
decay modes clearly support the molecular interpretation of the 7(3940) and 7(4140), while our estimates for the radiative 
decays provide a sensitive test for the underlying meson structure of the two 7 mesons discussed here. The alternative 
assignment of = 2++ is also tested, giving similar results for the strong decay widths. 
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INTRODUCTION 

With the discovery of the charmonium-likeX, Y and Z mesons the prospects were raised of possibly identifying a non- 
trivial structure of a meson resonance. The latest announcement [1] of the narrow state 7(4140) by the CDF Collabora- 
tion at Fermilab is just the continuation of a series of previously discovered charmonium-like resonances which are not 
easily interpreted as conventional quark-antiquark structures. Now the CDF Collaboration has evidence of a narrow 
near-threshold structure, termed the 7(4140) meson, in the J/ij/^ mass spectrum in exclusive Z?+ — > J/xj/^K^ decays 
with the mass otj'(4i4o) = 4143.0 ± 2. 9(stat) ± 1.2(syst) MeV and natural width rj.(4i4o) = 1 1.7+5 Q(stat) ± 3.7(syst) 
MeV [1]. As already stressed in [1], the new structure 7(4140), which decays to J/ just above the J/ y/^ threshold, 
shows a similar decay and production pattern to the previously discovered 7(3940) [2, 3], which is also produced in 
B-decays and was found in the J/ y/co decay channel near this respective threshold. The mass and width of the 7(3940) 
resonance are: my(394o) = 3943 ± 1 1 (stat) ± 13(syst) MeV, ry(394o) = 87 ± 22(stat) ± 26(syst) MeV (Belle Collabo- 
ration [2]) and mj.(394o) = 3914.6+|^(stat) ± 2.0(syst) MeV, TY(3g40) = 34+^2(stat) ± 5(syst) MeV (BABAR [3]). 

The problematic nature of these resonances with respect to their arrangement in the constituent quark model 
originates in the decay patterns. A conventional pure charmonium state decays dominantly to open charm modes 
while hidden charm decays would be strongly suppressed [4, 5] due the Okubo, Zweig and lizuka (OZI) rule. Under 
the assumption that the 7(4140) is a pure cc state a quantitative [6] estimate for the hidden charm 7(4140) J /W^ 
decay width results in the order of a few keV. This is in fact around several orders of magnitude smaller than current 
data seem to imply for the 7(4140) J / decay width. Similarly, the 7(3940) J / y/co decay width is larger than 
1 MeV [5], orders of magnitude away from the cc expectation. This could be signals for nonconventional structures 
of the 7(3940) and the 7(4140). Possible alternative interpretations involve structures such as hadronic molecules, 
tetraquark states or even hybrid configurations (for recent reviews see e.g. Refs. [4, 5]). 

The closeness of the D*D* and D*^D*^ thresholds implies a meson-meson bound state structure. But also quanti- 
tative analyses based on the possible dynamical generation of both resonances seem to support this interpretation. As a 
first follow-up to the CDF result it is suggested in [7] that both the 7(3940) and 7(4140) are hadronic molecules in this 
case in the context of a meson-exchange dynamics. These hadron bound states can have quantum numbers J^'" = 0++ 
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or 2++ whose constituents are the vector charm Z)*(Z)*) mesons: 



|y(3940)) = {\D*+D*-) + \D*°D*0)) , 
v2 

|y(4140)) = \D*+D*-). (1) 

Earlier results based on the pion-exchange mechanism already indicated that the D*D* system can form a bound 
state [8]. The dominant D*D* component in the 7(3940) is also confirmed in the coupled channel analysis of [9]. 
Binding in the D*D* channel can be induced by ?] and (j) meson exchange [7]. Further on recent QCD sum rule 
studies [10, 11] also favor a molecular structure of the 7(4140). 

In the present article we report on a first quantitative prediction for the decay rates of the observed modes 
7(3940) — > J/ij/CO and 7(4140) — t- J/y/(p assuming the hadronic molecule structures of Eq. (1) with quantum numbers 
/''^ = O''""'" (for further details see [12]). Results will be shown to be fully consistent with present experimental 
observations, strengthening the unusual hadronic molecule interpretation. Further predictions are given for the radiative 
two-photon decays of these states. Finally, we also consider the alternative = 2++ assigimient for the 7 states. 



THEORETICAL BACKGROUND 

In order to study the decay properties of hadronic bound objects we use a method based on effective Lagrangians 
describing the interaction of the particles involved. Other hadron molecules such as the ao(980), /o(980), Z)Jq(2317), 
Z)si(2460), and X(3872) have already been successfully described and studied with this method (see e.g. [13, 14]). 
The composite (molecular) structure of the 7(3940) and 7(4140) states is defined by the compositeness condition 
Z = [15, 16] which impUes that the renormalization constant of the hadron wave function is set equal to zero. 
Therefore the hadron exists solely as a bound state of its constituents. As a consequence decay processes of hadron 
molecules proceed in leading order via meson loops of its constituents. 

For the observed 7(3940) and 7(4140) states we adopt the convention that the spin and parity quantum numbers of 
both states are J^'^ = 0^'^. Presently the quantum numbers are not unambiguously determined yet in experiment 
except for C = +. For example, the 7(3940) is also discussed as a = l"*""*" charmonium candidate [5], but O"*""*" 
is not ruled out. Their masses are expressed in terms of the binding energy ey as ffJF(394o) = 2mo* — er(3940) 
"^7(4140) = 2mD* — ey(4i4o), where mo* = inu*+ = 2010.27 MeV and mo* = mjj*+ = 2112.3 MeV are the masses of 
the constituent mesons. Since the observed masses are relatively far from the corresponding thresholds we do not 
include isospin-breaking effects that is we suppose that charged and neutral nonstrange D* mesons have the same 
masses. Following Ref. [7] we consider the 7(3940) meson as a superposition of the molecular D*^D* and D*^D*^ 
configurations, while the 7(4140) is a bound state of D*'^ and D*~ mesons (see Eq. (1)). For the sake of simphcity we 
introduce the short notations 7(3940) = 7i and 7(4140) = 72. The coupling of the scalar molecular states Yj to their 
constituents is expressed by the phenomenological Lagrangians 

^y, = ^7i(x)/A^(y^){0*°''(x+|)D^^^(x-|)+Z)*+''(x+|)iP-(x-^)}, 

^Y, = gYjiix) j d'y^{y^)D*+^{x+^-)Dl-{x-^-), (2) 

where gy is the respective coupling constant.The interaction Lagrangian is expressed by the center of mass and relative 
coordinates x and y. The distribution of the constituents inside the molecular states 7, is expressed by the correlation 
function ^{y^). The Fourier transform of the correlation function appears as a form factor in our calculations. In the 
present analysis we have chosen the Gaussian form of ^{p\/J^) = exp(— p| /Ay), where pe is the Euchdean Jacobi 
momentum and Ky is a size parameter with a value of about 2 GeV - a typical scale for the masses of the constituents 
of the 7; states. 

The coupling constants gy are determined by the compositeness condition with Zy. = 1 — Ey, (my) = 0, where 
l,Y.{mY.) = dl,y.{p^)/dp^\p2^„2 is the derivative of the mass operator Ey,. generated by ^}^.(x). 

To determine the strong Y J / i/aV and two-photon 7 — )• yy decays we have to include the couphngs of D* (D*) 
mesons to vector mesons {J/ y/, (0, <j)) and to photons. The couplings of // i/A, o), <j) to vector D* (D*) mesons are taken 



from the HHChPT Lagrangian [17, 18]: 



■^D*D'Jxf, 

where AdB = AdB ~ BdA, i,j are flavor indices of the diagonal matrix y,y = diag{a)/\/2, fi)/\/2, 0} containing 
(O and mesons and D*^ ~ is the triplet of vector D* mesons containing light antiquarks u, d and 

s, respectively. The chiral couplings g^t^,, , ^o*d*v ^"'^ fn'o'v taken from [17, 18] The leading-order process 
relevant for the strong decays 7(3940) J/\lfCO and 7(4140) J /W^ is the diagram of Figure 1 involving the 
vector mesons D* or D* in the loop. Diagrammatically it already becomes obvious that the meson loop technique 
provides an opportunity to circumvent the OZI-suppression at lowest order in contrast to corresponding diagrams in 
the charmonium picture. 




FIGURE 1. Diagram describing the Y J/yV (V = CO, (j)) decays. 

The coupling of the charged D**(D**) mesons to photons is generated by minimal substitution in the free 
Lagrangian of these mesons. The corresponding electromagnetic Lagrangian is given by 

^em - eAa [g"' D*^- id^ D*+ - g^' D*' id'' D*+ + H.c) 

+ e^D*-D*+(A^'A^-g^^A"Aa). (4) 

which diagrammatically leads to the two relevant graphs displayed in Figures 2(a) and 2(b). Since the strong interaction 
Lagrangian of Eq. 2 is nonlocal it has to be modified in the following way in order to restore full gauge invariance. 
According to [19] each charged constituent meson field in Jiy is multiplied by the gauge field exponential 
^^(y) ~^ e^"''(^''^'^'//*(y), where I{x,y,P) — j!^,dZjj,A^{z)- This gives rise to the diagrams in Figures 2(c) and 2(d) 
which therefore are a consequence of the nonlocality of the present method. The contribution of these additional 
processes is of the order of a few percent when compared to the leading diagram of Figure 2(a). However, their 
inclusion is necessary to guarantee full gauge invariance. 




(e) (d) 

FIGURE 2. Diagrams contributing to the Y ^ yy decay. 



RESULTS 



In the present analysis we studied the strong J/xj/V (V = aJ,<^>) and radiative two-gamma decay properties of the 
y(3940) and 7(4140) mesons in the framework of the effective Lagrangian approach. Our results for the decay widths 
are contained in Table 1 . The errors of the results reflect the uncertainty of the experimental mass values of the Y 
states. For the masses of the Y mesons we use the values extracted by the BABAR |3] and the CDF [1] Collaborations. 

TABLE 1. Decay properties of y(3940) and y(4140) for 



Quantity 


Y(3940) 


Y(4140) 


r(y-)-7/v^y),Mev 


5.47 ± 0.34 


3.26 ± 0.21 


r(y -> 77), keV 


0.33 ± 0.01 


0.63 ± 0.01 



The predictions for the couplings gy of the Y states to their meson constituents are consistent with a trivial estimate 
using the Weinberg formula originally derived for the deuteron as based on the compositeness condition [15] with 

= \/ 327t m^J^ Sy^. This formula represents the leading term of an expansion in powers of the binding energy e. 
Note that this expression can be obtained in the local limit (i.e. the vertex function approaches the limit <E>()'^) — )• 5^(>')) 
and when the longitudinal part k^k^ /mj), of the constituent vector meson propagator is neglected. The numerical 
results for the couplings 8^1,940) = 9.16 GeV and ^7(4149) = 8-91 GeV are rather comparable with the nonlocal results 
^F(3940) = 14.08 GeV and ^^(4140) = 13.20 GeV. 

The predictions of r{Y(3940) J/y/co) = 5.47 MeV and r(F(4140) J/wP) = 3-26 MeV for the observed 
decay modes are sizable and fully consistent with the upper limits set by present data on the total widths. The result for 
r(F(3940) — )• // y/co) is also consistent with the lower Umit of about 1 MeV [5]. Values of a few MeV for these decay 
widths naturally arise in the hadronic molecule interpretation of the 7(3940) and F(4140), whereas in a conventional 
charmonium interpretation the J/y/V decays are strongly suppressed by the OZI rule [5]. In addition to the possibility 
of binding the D*D* and Z)*+D*" systems [7], present results on the J/ y/V decays give further strong support to the 
interpretation of the Y states as heavy hadron molecules. 

Further tests of the presented scenario concern the two-photon decay widths, which we predict to be of the order of 
IkeV 

TABLE 2. Decay properties of F(3940) and y(4140) for 
/^C = 2++. 



Quantity 


Y(3940) 


¥(4140) 


r{Y ^J/\lfV),MeV 


7.48 ±0.27 


4.41 ±0.16 


r(y 77), kev 


0.27 ±0.01 


0.50 ±0.01 



Finally we also test the = 2++ assignment for the Y states which is not yet ruled out experimentally. The 
coupling of the molecular tensor field Y^y to the meson constituents is set up as 

^y, = ^Yrix)Jd'y^{/){D;\.+ l)WOy{.-'-)+D;^x+'-)D;-{.-l)} 

= gy,Y,^\x)Jd'y<^{/)D*+{x+l)D*-{x-l) (5) 

Proceeding as outlined before we obtain the results in Table 2. Since the results for the strong J/y/ decays are quite 
similar to the O"*""*" case, a 2++ scenario cannot be ruled out and is also consistent within a molecular interpretation of 
the Y states. 

Quite recently the BELLE Collaboration [20] searched for the 7(4140) in the two-photon process yy — )• J/y/cj), 
however without evidence so far. In this framework they tested our prediction for the 7(4140) 77 decay width 
by experiment and obtained ryy(7(4140))^(7(4140) J/Y'P) < 40 eV [20] for = 0+ which results in a much 
smaller upper bound for the two-photon widths of about 0.2 keV. This finding is presently in conflict with a possible 
molecular interpretation of the 7(4140). 

The search for charmonium-like states in the process 77 J/yfca by the BELLE Collaboration [21] resulted in an 
enhancement with mass M = 3915 ± 3 ± 2 MeV denoted by X(3915). According to its mass and width it is a possible 



candidate for the 7(3940). The observed 2/ decay properties also support the D*D* bound state interpretation since 
the measured quantity, the product of the 27 width and the branching ratio, ryy(X(3915))^(X(3915) (OJ/^) = 
61±17±8eV (y^ = 0+) is of the same order of magnitude as the present results in Table 1. 

CONCLUSIONS 

The properties of the hidden-charm and two-photon decay modes of the charmonium-like mesons 7(3940) and 
7(4140) are discussed under the assumption that both are bound states of charm mesons. The method used in 
the present analysis constitutes a consistent field theoretical tool for hadronic bound states and fulfills full gauge 
invariance. In case of the hidden charm decays 7(3940) J/\f/(0 and 7(4140) J/yf'P the sizable results for the 
widths obtained in the effective Lagrangian approach are in good agreement with available data and clearly support 
the hadron molecule assignment. A cc scenario is disfavored since, in contrast to the hadron molecule, the estimated 
widths is negligible in this case. 

In addition the two-photon decays can provide a sensitive test for the meson structure as well. In case of the 27 
width of the 7(4140) available data predict a much smaller width. Note that if the present experimental upper limit 
for the 27 decay width is even lower it will be very hard to explain this state as a hadronic molecule. In contrast the 
7(3940) — > 77 decay width together with the sizable J/xj/O) branching ratio obtained within the D*D* bound state 
interpretation are comparable to experimental observations. 

However, more and particularly precise experimental information on these 7 states is necessary. A full interpretation 
of the 7(3940) and 7(4140) states requires an experimental determination of the quantum numbers, a consistent 
and hopefully converging study of binding mechanisms in the D*^^^D*^^-^ systems and finally a full understanding of the 

open charm decay modes, such as DD, DD*, DD*j, etc., which are also naturally fed in a charmonium picture. 
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